Cosmological constraints on Lorentz invariance violation in the neutrino sector 
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We derive the Boltzmann equation in the synchronous gauge for massive neutrinos with a deformed 
dispersion relation. Combining the 7-year WMAP data with lower-redshift measurements of the 
expansion rate, we give constraints on the deformation parameter and find that the deformation 
parameter is strong degenerate with the physical dark matter density rather than the neutrino mass. 
Our results show that there is no evidence for Lorentz invariant violation in the neutrino sector. 
The ongoing Planck experiment could provide improved constraints on the deformation parameter. 
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I. INTRODUCTION 

Neutrino oscillations imply that at least two of the 
three neutrino types have non-zero mass (see Ref. [1] 
and reference therein). Unfortunately, neutrino oscilla- 
tion experiments only provide two mass squared differ- 
ences, but not the overall mass scale. Cosmology provides 
a promising way to tackle this problem by the gravita- 
tional effect of massive neutrinos on the expansion his- 
tory near the epoch of matter-radiation equality and on 
the growth of large-scale structures at late time. Preci- 
sion measurements of the cosmic microwave background 
(CMB) anisotropics and the large-scale structure distri- 
bution of matter have allowed us to determine or con- 
strain the absolute neutrino mass scale [2-5] (see also 
Ref. [6] for a recent review). 

Indeed, neutrino oscillations can be explained by small 
Lorentz invariance violation even without introducing 
neutrino mass, as shown in [7, 8, 10-12]. Observed neu- 
trino oscillations may originate from a combination of 
effects involving neutrino masses and Lorentz invariance 
violation. The possibilities of Lorentz invariance viola- 
tion were explored in quantum gravity [13], loop quan- 
tum gravity [14], non-commutative field theory [15], and 
doubly special relativity theory [16]. 

It is well-known that Lorentz symmetry is a funda- 
mental feature of modern descriptions of the nature, in- 
cluding both the Einstein's general relativity and stan- 
dard model of particle physics. One can expect that 
the breaking of Lorentz symmetry may leave some im- 
prints in astrophysical observations such as the CMB 
anisotropics and large-scale structure of our universe. In 
this paper, we consider cosmological tests of Lorentz in- 
variance violence in the neutrino sector. Here we focus 
on the Coleman- Glashow model, in which the encrgy- 
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momentum relation is modified by a Lorentz-violating 
interaction in the framework of conventional quantum 
field theory [8]. We will build up a Lagrangian and de- 
rive the Boltzman equation in the synchronous gauge for 
the massive neutrinos with deformed dispersion relation, 
and then place constraints on Lorentz invariance viola- 
tion by combining the 7-year WMAP data [3] with the 
latest distance measurements from the Baryon Acoustic 
Oscillation (BAO) in the distribution of galaxies [17] and 
the Hubble constant (Hq) measurement [18]. 

This paper is organized as follows. In Sec. II we write 
down the Lagrangian for neutrinos with deformed disper- 
sion relation. In Sec. Ill we derive the Boltzmann equa- 
tion for neutrinos in the synchronous gauge. In Sec. IV 
we place constraints on the deformation parameters us- 
ing the CMB data in combination with measurements of 
the Hubble constant Ho and the BAO feature. Section V 
is devoted to conclusions. 



II. DEFORMED DISPERSION RELATION 

At a phenomenological level the deformed dispersion 
relation for massive neutrinos can be generally parame- 
terized by 

oo n 

E? =m *+f + Y l °nj^, (1) 
n— 1 

where E is the neutrino energy, m is the neutrino mass, 
P = {p l Pi) 1 ^ 2 the magnitude of the 3- momentum, a„'s are 
dimensionlcss coefficients and M denotes the energy scale 
corresponding to Lorentz symmetry violation (which is 
typically taken to be the Planck mass). Such a deformed 
dispersion relation implies that there are departures from 
Lorentz invariance in the neutrino sector if a n ^ 0. The 
n = 1 term would produce huge effects at low energy and 
has been strongly constrained. The p n term with n > 3 
is suppressed by 1/M n ~ 2 . In the present work we will 
therefore consider the case of n = 2, 

E 2 = m 2 + p 2 + £p 2 , (2) 
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where £ is the deformation parameter characterizing the 
size of Lorentz symmetry violation. Such a deformed 
dispersion relation was constructed by Coleman and 
Glashow in the framework of conventional quantum field 
theory [8]. 

Here we point out that the dispersion relation given in 
(1) and (2) is not very general. It neglects oscillations, 
possible species dependence, anisotropics associated with 
violations of rotation symmetry, and CPT violation. As 
shown recently by Kostelecky and Mewes, all of these are 
possible [9]. For example, it is known that taking odd 
values of n in Eq. (1) corresponds to CPT violation and 
gives rise to a sign difference for neutrinos and antineutri- 
nos [9] . We have to emphasize that the model considered 
in this paper is one of many possible Lorentz- violating 
theories. 

The scalar perturbations of the Fricdmann-Lcmaitre- 
Robertson- Walker metric in the synchronous gauge can 
be written as 



ds 2 = a 2 (r) [-dr 2 + {5 l} + h l3 )dx l dx j ] , 



(3) 



where a(r) is the scale factor, r is the conformal time, 
and the scalar mode of hij is represented by two functions 
h and rj which are defined by 



h i:j (x,T) = / d 3 ke iks 



in Fourier space, where k = kk. The action for a neutrino 
with the dispersion relation (2) can be written in the first 
order formalism by 

S = Jdr [±"P M - A (m 2 + g^P^ + ZgVPiPj)] , (5) 

where A is a Lagrange multiplier which enforces the mass- 
shell condition (2) and 
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The above action implies that the Hamiltonian van- 
ishes. Performing the Legendre transformation to the 
Lagrangian formalism and eliminating A from the action 
(5), we obtain the Lagrangian for neutrino as follows 



-ma* / 1 
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hijU l u J 



2(l + £-u 2 ) 



(8) 



where u l = dx 1 jdr is the coordinate velocity and u 2 = 
5ijU l u^ . Following [19] from the Lagrangian (8), the 
equations of motion in the synchronous gauge become 



dx l 

~dT 

dq 



qhijn l n J , 



[l + 0(h)\, (9) 
(10) 



where qj = qrij is the comoving 3-momentum written 
in terms of its magnitude and direction with n l ni = 
dijifni = 1, and e = y m 2 a 2 + (1 + £)q 2 is the comoving 
energy of neutrinos. 



III. BOLTZMANN EQUATION 

The number density n Ul energy density p v and pressure 
P v for massive neutrinos with (2) are given by 
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fo(q) , (11) 

^fo(q) , (12) 
£)<7 2 

— Jo(<z) ■ (13) 



Here the zeroth-ordcr phase space distribution is well ap- 
proximated by the relativistic Fermi-Dirac distribution 
fo(q) = .9s [1 + cxp(v"l + tq/To)]- 1 , where T is the neu- 
trino temperature today and g s — 2 the number of spin 
degrees of freedom. If to ^ To, the total mass of neutri- 
nos m = 94(1 + £)3/ 2 ft„/i 2 cV. 

The perturbed energy density, pressure, energy flux, 
r£\ and shear stress in the Fourier space k are respectively 
given by 



= A / tS^/o(<7)*o, (14) 
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where the perturbations *f?i satisfy the following Boltz- 
mann equation 



2 f d 3 » 
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in the synchronous gauge. Such a Boltzmann hierar- 
chy is effectively truncated by adopting the following 
scheme [20] 
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The initial conditions for the perturbation vp; and r\ on 
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super- horizon scales (fcr <C 1) are given by 

W = rO, 



4 " ding ' 
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(21) 
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r, = 2C- 5 + 9VTT ^ 5 ^C( fc rf, (23) 
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where 



K = --C(kT) 



(24) 
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C(fcV). (26) 



Here, C is a dimensionless constant determined by the 
amplitude of the fluctuations from inflation and R u = 
Pv/(pv + Pi) during radiation domination. 

In order to compute the theoretical CMB power spec- 
trum, we modified the Boltzmann CAMB code in [27] 
to appropriately incorporate the Lorentz-violating term 
in the neutrino sector. Actually the Lorentz-violating 
term affects not only the evolution of the cosmological 
background but also the behavior of the neutrino pertur- 
bations. From Eqs. (11)-(13), we see that increasing £ 
decreases the number density, energy density and pres- 
sure of neutrinos, and thereby increases the redshift when 
the matter density equals the radiation density and re- 
duces the expansion rate prior to and during the epoch of 
photon-baryon decoupling. This leads to reduced heights 
of the first and second peaks of the CMB while a nearly 
constant increase in the acoustic oscillation amplitudes 
at I > 600. On the other hand, the coefficient of the sec- 
ond term in the Boltzmann equation (18) plays an active 
role, which alters the shape of the CMB power spectrum 
caused by changing neutrino propagation. Decreasing £ 
increases fluctuation power both at I < 10 and at I > 100. 
Moreover, the CMB is more sensitive to negative values 
of £ than positive ones. These two effects can be distin- 
guished from a change in the total mass of neutrinos or in 
the effective number of extra relativistic species [21-26], 
as shown in Figure 1. 



IV. COSMOLOGICAL CONSTRAINTS 

In our analysis we use a modified version of the pub- 
licly available CosmoMC package to explore the param- 
eter space by means of Monte Carlo Markov chains tech- 
nique [28]. We consider a flat ACDM models plus three 
Lorentz-violating neutrino species, described by a set of 
cosmological parameters 

{fl b h 2 , il c h 2 , 9 S , r, n s ,A s ,T,m, £}, 




FIG. 1: Theoretical angular power spectrum of the CMB for 
£ = -0.1, 0, 0.1. Here £ m = 3 x 0.3 eV is fixed. 



where h is the dimensionless Hubble parameter such that 
Hq = 100/i kms -1 Mpc -1 , Vt b h 2 and fl c h 2 are the physi- 
cal baryon and dark matter densities relative to the crit- 
ical density, S is the ratio of the sound horizon to the 
angular diameter distance at the photon decoupling, r is 
the reionization optical depth, n s and A s are the spectral 
index and amplitude of the primordial curvature pertur- 
bations at the pivot scale ko = 0.002 Mpc -1 , is 
the total mass of neutrinos assuming that the three neu- 
trino masses are approximately degenerate, and £ is the 
deformation parameter. 

We use the seven-year WMAP (WMAP7) data with 
the likelihood code supplied by the WMAP team. We 
consider the Sunyaev-Zel'dovich (SZ) effect, in which 
CMB photons scatter off hot electrons in galaxy clusters. 
Given a SZ template, the effect is described by a SZ tem- 
plate amplitude Agz as in the WMAP paper [3] . We also 
impose Gaussian priors on the Hubble constant, Hq = 
74.2 ± 3.6 kms -1 Mpc -1 , measured from the magnitudc- 
redshift relation of low-z type la supernovae [18], and on 
the distance ratios, r s /D v (z = 0.2) = 0.1905 ± 0.0061 
and r s /D v (z = 0.35) = 0.1097 ± 0.0036, measured from 
BAO in the distribution of galaxies [17]. Here r s is the 
comoving sound horizon size at the baryon drag epoch 
and D\r is the effective distance measure for angular di- 
ameter distance. 

In Table I we summarize our results. With the data 
of WMAP7+i?o+BAO, the deformation parameter is es- 
timated to be £ = —0.077 ± 0.089, which implies a null 
detection of Lorentz invariant violence within the error 
limits. Such large uncertainties in £ mainly come from 
the strong correlation between the deformation param- 
eter and the physical dark matter density as shown in 
Figure 2. In this case the uncertainties in il c h 2 are about 
three times larger than those derived in the standard 
ACDM model [3]. Since the deformation parameter is 
nearly uncorrelated with the total mass of neutrinos in 
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parameter WMAP+Ho+BAO WMAP+iJp+BAO Planck+iJp+BAO 

W0Q b h 2 2.312 ±0.078 2.314 ± 0.081 2.260 ± 0.023 

Q, c h 2 0.1200 ±0.0103 0.1206 ±0.0104 0.1133 ± 0.0033 

Q A 0.715 ± 0.020 0.723 ±0.017 0.725 ± 0.005 

H 70.9 ± 2.4 km/s/Mpc 72.1 ± 2.2 km/s/Mpc 70.3 ± 1.0 km/s/Mpc 

r 0.089 ±0.014 0.088 ± 0.014 0.088 ± 0.004 

£m < 0.51 eV (fixed) (fixed) 

£ -0.077 ±0.089 -0.073 ± 0.089 -0.005 ± 0.037 



TABLE I: Mean values and marginalized 68% confidence level for the deformation parameter and other cosmological parameters. 
For the total mass of neutrinos, the 95% upper limit is given. 



Figure 2, our constraints on the deformation parameter 
are not significantly changed if three neutrinos are mass- 
less, as we can see from Table I. Compared to parti- 
cle physics experiments, cosmological observations yield 
much weaker constraints on the Lorentz-violation param- 
eter in the neutrino sector. As listed in Table XIV of 
Ref. [29], previous constraints range from parts in 10 5 
to parts in 10 15 from time-of-flight measurements and 
various threshold analyses. 

Moreover, we present the constraints on the deforma- 
tion parameters from the ongoing Planck experiment [30] 
in Table I. Following the MCMC method described in 
Ref. [31], we generate synthetic data for the Planck ex- 
periment and then perform a systematic analysis on the 
simulated data. As we can see in Table I, the Planck data 
plus measurements of the Hubble constant and the an- 
gular diameter distance will reduce the uncertainties in £ 
by a factor of 2.4. Therefore, Planck CMB measurements 
allow us to detect the signature of Lorentz invariant vio- 
lence at 2er confidence level if |£| > 0.074. 

V. CONCLUSIONS 

We studied the cosmological consequences of Lorentz- 
violating neutrinos. We obtained the generalized La- 



grangian for a neutrino with the Coleman-Glashow type 
dispersion relation, and then derived its Boltzmann equa- 
tion and initial conditions in the linear theory of cosmo- 
logical perturbations. Using the 7-ycar WMAP data in 
combination with measurements of the Hubble constant 
and the BAO feature, we found no evidence for Lorentz 
invariant violence in the neutrino sector. The ongoing 
Planck experiment is expected to able to give a more 
stringent constraints on the deformation parameter. 
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